Motion of vortex lines in quantum mechanics* 
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Exact analytic solutions of the time dependent Schrodinger equation are produced that exhibit 
a variety of vortex structures. The qualitative analysis of the motion of vortex lines is presented 
and various types of vortex behavior are identified. Vortex creation and annihilation and vortex 
interactions are illustrated in the special cases of the free motion, the motion in the harmonic 
potential, and in the constant magnetic field. Similar analysis of the vortex motions is carried out 
also for a relativistic wave equation. 



I. INTRODUCTION 



Vortices have been a source of fascination since time 
immemorial. Empedocles, Aristotle pL and Descartes 
Q tried to explain the formation of the Earth, its grav- 
ity, and the dynamics of the whole solar system as due 
to primordial cosmic vortices. Lord Kelvin (^] envisaged 
atoms as vortex rings. Modern science of hydrodynamic 
vortices began with the fundamental work of Hclmholtz 
j| on the theory of vorticity. The present knowledge of 
vortex dynamics in fluid mechanics is best summarized 
in a recent monograph by Saffman Q . 

In our paper we shall be concerned with vortices not in 
ordinary fluids but in a peculiar fluid — the probability 
fluid — that can be associated with the wave function of a 
quantum system. Vortices of that type have been studied 
in the past in connection with superfluidity, superconduc- 
tivity, and recently in connection with Bose-Einstein con- 
densation |^-^] . The main tool in the study of these phe- 
nomena is the Ginzburg-Landau (Gross-Pitaevskii) equa- 
tion. The nonlinear character of this equation greatly 
complicates the analysis of time evolution of the vor- 
tices, especially in three dimensions. There are no exact 
analytic solutions of this equation with nontrivial vor- 
tex structure. One must resort either to approximations 
which can reveal only some features of the vortex motion 
(as, for example, in the early work of Fetter 10 or to 
extensive numerical calculations (cf., for example, [^2| ). 
In the present paper we use the standard time-dependent 
Schrodinger equation to investigate the motion of vortex 
lines embedded in the probability fluid of a quantum par- 
ticle. The linearity of this equation enables us to obtain 
a variety of exact, time-dependent analytic solutions for 
the wave functions that describe the probability flow with 
one, two or more moving vortex lines. We are also able 
to study the motion of vortex lines in the presence of ex- 



ternal forces. Our analysis of many special cases makes 
it possible to discover certain general phenomena such 
as the vortex switchover and the vortex creation and an- 
nihilation. The important conclusion of our analysis is 
that quantum effects play a crucial role in those phenom- 
ena. We believe that our analysis captures the essential 
features of quantum vortex dynamics that are common 
to many different situations. In particular, our analysis 
also extends to the relativistic domain. 

We have to admit that vortex lines associated with 
wave functions are rather elusive objects. They might be 
connected with physical reality in two ways. First, the 
vortex lines can appear only at those places where the 
wave function is equal to zero; the vanishing of the prob- 
ability density definitely has observable consequences. 
Second, the wave function of a macroscopically populated 
state of many bosons acquires a direct physical interpre- 
tation; vortex lines in this case m ay bec ome directly ob- 
servable as, for example, in BEC 
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IL VORTEX LINES IN THE HYDRODYNAMIC 
FORMULATION OF WAVE MECHANICS 

In quantum theory vortex lines arise in the hydrody- 
namic formulation of wave mechanics. This formulation 
originally due to Madelung |ll| (see also a review arti- 
cle ]l6| ) employs the hydrodynamic variables to describe 
the flow of quantum mechanical probability for a single 
particle. The probability density p(r, t) and the velocity 
field v(r, t) are connected with the wave function ^(r, t) 
through the formulas 



tp(r,t) = R(r,t)exp[iS(r,t)/h] 



p(r,t) = \^(r,t)\ 2 = R 2 (r,t) ) 



(1) 



(2) 
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_ 1 Re[ip*(T,t)(-iKV - eA)V>(r,t)] 
V(r '^ ~ m |V(r,t)| 2 

= -(V5(r,t)-eA(r,t)), (3) 
m 

where m is the mass of the particle and A(r, f) is the 
electromagnetic vector potential. The gradient term VS 
in (||) does not contribute to vorticity. Thus, the bulk 
vorticity V x v(r, t) in the probability fluid is completely 
determined by the magnetic field except at the points 
where the phase S is singular. This may occur only at 
the points where the wave function vanishes. The van- 
ishing of the complex wave function, in turn, gives two 
conditions: the vanishing of two real functions. Each 
condition defines a two-dimensional surface. The inter- 
section of these two surfaces defines a line. Thus, in 
addition to the given distribution of vorticity in the bulk 
of the probability fluid, that is uniquely determined by 
the magnetic field, we may also have isolated vortex lines 
embedded in the fluid. Along these lines the vorticity has 
a two-dimensional Dirac delta- function singularity in the 
plane perpendicular to each vortex line, ft is the dynam- 
ics of these isolated vortex lines that will be the subject 
of our study. In this context we shall identify vortex lines 
with line vortices. 

To ensure that the wave function is single valued, the 
strength of every vortex as measured by the circulation 
r along any closed contour C encircling the vortex line, 

r = £dl-v(r,t), (4) 

must be quantized: 

r= — n, n = 0,±l,±2,.... (5) 

TO 

As one approaches the vortex line the velocity v(r, t) 
must tend to infinity in order to satisfy the quantization 
condition (||). 

The circulation is conserved during the time evolution 
jl7],[l8| (in full correspondence with the Helmholtz-Kelvin 
law for ordinary nonviscous fluids p9|). The creation of 
vortex lines may proceed through two different scenarios: 
they may either appear in the form of a closed vortex line 
that springs from a point (Fig. |l|) or they can be created 
in pairs of opposite circulation (Fig. ^ . 

The hydrodynamic formulation can also be introduced 
in wave mechanics of a sp inning particle governed by the 
Pauli equation 

Despite all the formal similarities between the hydro- 
dynamics of classical fluids and the hydrodynamic for- 
mulation of quantum mechanics, there is a fundamental 
difference between the properties of line vortices in these 
two cases. In classical hydrodynamics, as pointed out by 
Saffman j5|, p. 36, "the line vortex is a useful kinematic 
concept, but it does not exist as a sensible dynamical 
limit" . In quantum mechanics the line vortices occur in 
exact solutions of the Schrodinger equation without any 



limiting procedure. The fact that the velocity becomes 
infinite as one approaches the line vortex is acceptable, 
because it is only the probability and not any real matter 
that flows with that velocity. 



III. THE ANATOMY OF A VORTEX LINE 

The commonly known examples of vortex lines are 
those embedded in the wave functions that have a fixed 
value hm z of the z-component of angular momentum 
(eg., the spherical harmonics). These vortex lines run 
from — oo to oo along the z-axis and have the circula- 
tion Hm z /m. Usually they are stationary, due to angular 
momentum conservation. In addition to straight vortex 
lines, we can also easily produce stationary vortex rings 
by superposing wave functions with the same energy but 
different angular momenta. The simplest example is the 
superposition (with a ir/2 phase difference) of the elec- 
tronic ^200 and ip2io states in hydrogen. The vortex ring 
is formed in z = plane and has the radius of two Bohr 
radii. 

In our study we shall discuss a broader class of vortices 
that, in general, do not exhibit cylindrical symmetry even 
when they form a straight line. Let us consider an arbi- 
trary vortex line £(s) parametrized by the length s along 
the line. A typical quantum vortex is associated with the 
line of first order zeros of the wave function. Higher or- 
der zeros will, in general, lead to vortices that carry more 
that one unit of vorticity. The local properties of such a 
vortex at a point £(s) are determined by the behavior of 
the wave function near this point 

«(r -€(*))• VV(€(*)), (6) 

where we have kept only the lowest terms in the Taylor 
expansion. Since quantum vortices have zero thickness, 
this approximation is always valid sufficiently close to the 
vortex line. At each point on the vortex line let us de- 
fine a complex vector w(s) — the gradient of the wave 
function 

w(«) = VV(r)| p=€( , ) . (7) 

We shall show that this vector describes the basic prop- 
erties of the vortex. Let us note that since the derivative 
of the wave function along the vortex line with respect to 
the length parameter s is zero, the real and the imaginary 
parts of the vector w(s) lie in the plane perpendicular to 
the tangent vector t(s) = d£(s)/ds, 

— = t(s) -w(s) = 0. (8) 

The general complex vector w has 4 real components, 
but only 2 parameters contain the information about the 
vortex, since the multiplication of the wave function by a 
complex number does not change, of course, the velocity 
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field. Indeed, the velocity field v(x,y,z) calculated for 
the approximate wave function (p), 



h I w w 

v{x,y,z) = — : — - - 



2mi \wr w*r 

H r x (w x w*) 
2mi |w-r| 2 



(9) 



is a homogeneous function of w. Therefore, v does not 
change when w is multiplied by any complex number. 
For simplicity we have chosen the origin of the coordinate 
frame at the point £(s). Only in the special case, when 
the real and imaginary part of w are mutually orthogonal 
and of equal length, the velocity lines follow circles in the 
plane orthogonal to the vortex line. In all other cases the 
vortex is "squeezed" ; the velocity lines follow ellipses. In 
the degenerate case, when the real part and the imagi- 
nary part of w are parallel, the ellipse is squeezed into a 
line and the velocity disappears. The vector w will, in 
general, vary as we move along the vortex line. However, 
one may check by a direct calculation that the circulation 
([|) calculated for the velocity field (^) does not depend 
on w and is equal to ±27r?i/m (the sign depends on the 
orientation of the contour C) . Higher values of the cir- 
culation are obtained when the wave function has zeros 
of the higher order but the typical case, of course, is the 
first order zero. 

The above analysis shows that the appearance of vor- 
tex lines in wave mechanics does not require that some 
special conditions be met. On the contrary, a generic 
wave function has vortex lines. The vanishing of the real 
and the imaginary part of the wave function defines a 
line. The Taylor expansion of the wave function around 
the point lying on this line defines the complex vector w 
that determines the local structure of the vortex. 

The vector w not only determines the velocity field 
but also plays a crucial role in determining the motion 
of the vortex line itself. The velocity u = d£(s,t)/dt of 
a point £(s,t) lying on the vortex line can be obtained 
from the condition 



dip(£(s, t),t)/dt = u-w + dtp/dt = 0. 



(10) 



The vector u is defined only up to a vector parallel to the 
vortex line and it can be determined by solving Eq.(10): 



w x w 



u 



dip* 



dip 
— i 

dt 



Using the Schrodinger equation, 



we obtain 



^(r,t), 



2mi w x w* 



(wAf + w*Ai/>) . 



(11) 



(12) 



(13) 



Let us notice that the potential term V does not appear 
in the formula for the velocity u since the wave function 
vanishes on the vortex line. The motion of the vortex 
line at a given point is completely determined by the lo- 
cal properties of the wave function: its gradient and its 
Laplacian at this point. There is no direct relation be- 
tween the velocity u of the vortex line and the velocity 
v of the probability fluid. 



IV. GENERAL METHOD OF GENERATING 
VORTEX LINES 

Let us choose the initial condition for the solution of 
the time-dependent Schrodinger equation in the form 



0k(r) = exp(ik-r)0 o (r). 



(14) 



By differentiating <p^ with respect to the components of 
the vector k any number of times we obtain new wave 
functions. Each differentiation brings down a component 
of the position vector and in this manner we may gener- 
ate an arbitrary complex polynomial that multiplies the 
initial wave function. Carrying out the differentiations, 
adding the results with appropriate complex coefficients, 
and setting k = at the end, we arrive at the expression 
for the initial wave function of the form 



(W fl (r) + iW>(r))to(r), 



(15) 



where Wr and Wj are real polynomials in the three vari- 
ables x, y, and z. 

The conditions Wr(t) — and W/(r) = define two 
two-dimensional surfaces in the three-dimensional space. 
When these surfaces intersect, they define a line. This 
will be a vortex line if the circulation around it is dif- 
ferent from zero, which is the generic case. This proce- 
dure leads to wave functions that have initially vortex 
lines of almost arbitrary shape and topology embedded 
in them. In other words, the wave function 0k serves as 
the generating function for vortex lines. In practice we 
are facing technical limitations; for higher order polyno- 
mials the analysis becomes very cumbersome. However, 
if the polynomial Wr + iWi is constructed as a product 
of some simple polynomials in the form 

W H (x,y,z) + iWi(x,y,z) = {W' R (x,y,z) + iW' I {x,y,z)) 
x (W£(x,y,z)+iW?(x,y,z))..., (16) 

then the initial wave function is easily analyzed: it has 
several vortex lines, each corresponding to one of the fac- 
tors in (flit). Most of our discussion will be restricted 
to second degree polynomials which already lead to a 
plethora of different shapes and behaviors. 

The time evolution of vortex lines that were embedded 
in the initial wave function is determined, of course, by 
the Schrodinger equation. The advantage of using the 
generating function (|lj) is that once we find the time- 
dependent solution i/)k(r,t) satisfying the initial condi- 
tion 
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V> k (M = 0) =0 k (r), 



(17) 



we may use it to generate all solutions that at t — have 
the form (^) . Since the vector k does not appear in the 
Schrodinger equation, the differentiations with respect 
to its components do not interfere with the time evolu- 
tion: k-derivatives of the time-dependent wave function 
are also solutions of the Schrodinger equation. Carrying 
out the differentiations with respect to the components 
of k on the time dependent function ipk(r,i), we shall 
obtain the time evolution of vortex lines from the time 
dependence of the polynomials Wn(r,t) and Wi(r,t) in- 
troduced in the formula (|l5|). 

Obviously, wave functions are not restricted to those 
that have only a finite number of vortex lines described 
by a polynomial prefactor (|l5|). For example, replacing 
the sum of polynomials by the sum of trigonometric func- 
tions cos(x/a) +icos(y/b) leads to an infinite "forest" of 
vortex lines pointing in the ^-direction. 



V. VORTEX LINES FOR A FREELY MOVING 
PARTICLE 

In this Section we shall describe the time evolution of 
vortex lines for a freely moving particle. The simplest so- 
lution of the corresponding time-dependent Schrodinger 
equation, that has at t = the form (p^), corresponds to 
Mr) = 1, 



ipk(r, t) = exp(ik-r) exp(— ihk 2 t/2m). 



(18) 



Such a plane-wave solution is not square-integrable but 
that can easily be corrected by multiplying ( |l8| ) by a 
slowly varying envelope without affecting significantly 
the dynamics of vortex lines. We shall return to this 
problem at the end of this Section. 

Let us begin with the simplest solution that initially 
has one rectilinear vortex described by the wave function 
(xcosx + iysinx) exp(ik-r). The parameter x measures 
the degree of squeezing of the velocity field. Differentiat- 
ing the wave function ( |Tg| ) with respect to k x and k y we 
obtain the wave function with a vortex line. Since the 



parameters k x ,k y , and k z in (18) can be interpreted as 



the components of the particle momentum, hk — p, it is 
instructive not to set k to zero but keep it finite after the 
differentiation. This will result in a free motion of the 
whole vortex line with the velocity v = Kk/m: 



{x(t) cos x + iy{t) sin X ) V>k(r, t) , 



(19) 



where x(t) — x — v x t and y(t) = y — v y t. Throughout 
this paper x(t),y(t) and z(t) will always denote a free 
classical motion with constant velocity. Eq. ([ll]) means 
that the vortex line riding on the plane wave is moving 
uniformly, without changing its shape, with the speed 
determined by the particle momentum. This was to be 
expected as a consequence of the Gallilean invariance of 



the theory. The same free motion of a single vortex line 
was found in the nonlinear Schrodinger equation fLOf . 

Our next example of a vortex is a ring of radius R. It 
can be obtained, for example, by cutting a cylinder by a 
plane: 

W R (x,y,z) = x 2 +y 2 -R 2 , Wi(x,y, z) = az. (20) 

The ratio of the two dimensional parameters R and a de- 
termines the degree of squeezing of the vortex, tan x — 
a/2R. The motion of this vortex is described by the 
time-dependent wave function 



V> C M) = 



:(t) 2 + y{t) 2 -R 2 + ia[z(t) + — ) 

ma J 



x Vk(r,t). 



(21) 



The quantum correction (the term 2ht/ma) to the clas- 
sical vortex motion that appears here amounts only to a 
change of the vortex velocity in the direction of the cylin- 
der axis. In this example, the vortex ring in the form of a 
circle has been produced by intersecting a cylinder and a 
plane. However, there are many other possibilities, even 
with the use of only simple surfaces. It is easy to ob- 
tain at t = the same shape of the vortex line choosing 
different wave functions. The motion of these vortices 
may change dramatically, when one pair of surfaces is 
replaced by a different one. The circular vortex formed 
by intersecting a cylinder with a plane, as shown above, 
moves uniformly along the cylinder axis. On the other 
hand, the same vortex produced at t — by a sphere 
and a plane will contract to a point and disappear. The 
appropriate wave function has the form 



* s (r,0 



ma 



(22) 



The circle disappears, i.e., the vortex is annihilated, at 
the time t a — maR/Zh. Extending our analysis to nega- 
tive values oft, we see that this vortex was created at the 
time t c = —maR/2>h (Fig. |l|). A circular vortex obtained 
by intersecting a paraboloid with a plane has again a dif- 
ferent motion: after it is created, it keeps expanding and 
moving along the z-axis to infinity. Other motions of a 
circular vortex line can also be obtained by choosing a 
second quadratic surface, instead of a plane. In all these 
cases the motion of the vortex ring has two components. 
One is the free motion with the velocity v and the other 
is a quantum correction that may lead to the vortex an- 
nihilation. We would like to stress that the creation and 
annihilation of vortex rings are of a purely quantum ori- 
gin. These phenomena are not found in the studies of 
motion of the vortex ring based on the semiclassical ap- 
proximation |jniyri| . 

The creation and annihilation of vortex rings at 
a point, shown in Fig. [j], does not contradict the 
Helmholtz-Kelvin theorem on the conservation of circu- 
lation. This theorem, obviously, holds only inside the 
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fluid but the vortex creation and the vortex annihilation 
occur at the points where the fluid density vanishes. 

The mutual interaction between vortex lines can be 
seen in the motion of two vortex lines. The initial wave 
function containing two arbitrary rectilinear vortices has 
the product form 

0k (r) = (wi-(r - ri))(w 2 -(r - r 2 )) exp(ik-r). (23) 

The corresponding time-dependent wave function, ob- 
tained by our method is 

iht 

^(r,t) = [(wi-(r(£) - ri))(w 2 -(r(t) - r 2 )) H -w r w 2 ] 



x 1/Jk(r,t). 



(24) 



Thus, the two vortices become, in general, entangled by 
the quantum correction term: the wave function is not a 
product of factors describing two separate vortices. The 
coupling between the two vortices is due to the inter- 
action between the vortices and it vanishes only when 
wi-w 2 = 0. The time evolution exhibiting the phe- 
nomenon of a switchover of the two vortex lines is shown 
in Fig. [| To simplify the calculations we assumed that 
both vortices are nonsqueezed. The coordinate system is 
chosen in such a way that at t — the vortex lines lie in 
the z — ±a planes with the angle between them equal to 
2p. With this choice, the time-dependent wave function 
can be written in the form 

tp(r, t) = [WiW 2 - (2iht/m) sin 2 <p] tp k (r, t), (25) 

where the two polynomials W\ and W 2 , 

Wi(r, t) = x(t) cos ip + y(t) sin tp + i(z(t) + a), (26) 
W 2 (r, t) = x(t) cos tp - y(t) sin ip + i(z(t) - a), (27) 

describe the motion of both vortex lines separately, 
whereas the last term in the square bracket is respon- 
sible for their interaction. Only when ip — the vor- 
tex lines remain unchanged during the time evolution. 
This behavior is different from what can be found in 
non-compressible fluids p9| . The quantum interaction 
between the vortex lines is the strongest when the lines 
are antiparallel (i.e. the two vortices have opposite circu- 
lation, ip — 7r/2). In this case the vortex lines retain their 
direction but the distance between them decreases and at 
time t a = ma 2 /h the vortex lines collide and annihilate. 
The same wave function describes at t c — —ma 2 /fi the 
creation of two rectilinear vortices of opposite circulation. 
These processes are shown in Fig. ||. 

We shall now replace the plane wave solution of 
the Schrodinger equation, corresponding to 0o( r ) = 
1, by a normalizable Gaussian wave packet 4>o( v ) = 
cxp(— r 2 /2l 2 ). This will give us square integrable solu- 
tions with vortex lines. The generating function for such 
solutions has the form: 



exp(-k 2 P/2) 
{l+iht/ml 2 )i 



■ exp 



-(r-ikl 2 



2l 2 (l + iht/ml 2 ) 



(28) 



We shall illustrate the influence of the shape of the Gaus- 
sian envelope on the motion of vortices in the simplest 
case of one rectilinear nonsqueezed vortex that is paral- 
lel to the z-axis and passes through the point (xq , 0) at 
t = 0. The motion of such a line is composed of the free 
motion with the velocity v of the center of the Gaussian 
envelope and the motion induced by the shape of the en- 
velope. The time-dependent wave function in this case 



x(t) — xq + i(y(t) — xoht/ml 2 ) 



1 



(29) 



iht /ml 2 

The last term in the numerator represents the quantum 
correction to the free vortex motion and it vanishes in the 
limit, when the width of the Gaussian envelope I tends to 
infinity. The speed of this additional motion depends on 
the smallest distance of the vortex line from the center 
of the Gaussian but does not depend on the motion of its 
center. 



VI. VORTEX LINES FOR A PARTICLE IN A 
UNIFORM MAGNETIC FIELD 

In order to determine the motion of vortex lines in a 
magnetic field we need the counterpart of the generating 
function ( [l8|) in the presence of the field. We shall use 
the Schrodinger equation in the symmetric gauge 

*2 



id t ip(r,t) 
iheB 



{xd y 



2m 

yd x ) + ^(x 2 +y 2 )]^(r,t), (30) 



e 2 B 2 



x exp[?i(e 
x exp[i(e~ Mc * 

-itu c t 



m 8m 
where the z-axis is taken along the magnetic field. The 
solution of this equation, that at the time t = has the 
form exp(ik-r) exp(— eB(x 2 + y 2 )/4h) (a plane wave mul- 
tiplied by the ground state Gaussian in the xy-plane), 
is 

^ (r, t) = exp{-eB(x 2 + y 2 )/4h) exp[-iu c t/2] 

UA -l)(k 2 x + kl)/(2eB)] 

l)(xk x +yk y )/2] 

x exp[(e-^ J ~ l)(xk y ~ yk x )/2] 

x exp[izk z - ihk 2 z /(2eB)}, (31) 

where uj c = eB/m is the cyclotron frequency. Since the 
differentiation with respect to the components of k brings 
down linear terms in the space variables, a single vortex 
along a straight line will preserve its form. All that can 
happen is the motion of this straight line. For example, 
a line vortex that at t — lies in the y — a plane at an 
angle ip to the field direction will move in time according 
to the following parametric representation 

2a + xsm(uj c t) (1 + simp) 



y{x) 



z(x) = — 



1 — sin ip + (1 + sin ip) cos(uj c t) ' 
2x tan tp + a sin(w c i) (sec ip + tan ip) 



1 — sin 99 + (1 + simp) cos(uj c t) 



(32) 
(33) 
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In Fig. ^ we show the motion of this line during one pe- 
riod of the cyclotron motion. 



VII. VORTEX LINES FOR A PARTICLE IN A 
TRAP 

The study of vortex dynamics in a trap has a special 
significance because it may, perhaps, throw some light 
(as the first, linear approximation) on the behavior of 
the atoms that form the Bose-Einstein condensate. For 
simplicity, we shall consider a trap in the form of a spher- 
ically symmetric harmonic oscillator with the frequency 
u> and choose as </>o( r ) m the formula ( |l4| ) the ground 
state wave function exp(— mojT 2 / 2fi) . The time depen- 
dent wave function that serves as a generator of vortex 
lines in this case is 

ip k (r, t) = exp(-3zwi/ 2 ) ex P(~~^~ r ) 

x exp ( ze- lwt (k-r - ^— sm(ut)) \ . (34) 
\ 2muo J 

As an example we shall choose the vortex line that at 
t = has the form of a circular ring of radius R pass- 
ing through the center. The wave function describing the 
time evolution of such a vortex is obtained by evaluating 
the appropriate derivatives of (^i|). Setting k = after 
the differentiation, we obtain 



[e *—(.»- 



-IlLOt / 2 



y 



— ) + — 



- e~ iut R{2x - iz)} 
xCo(r.t)- (35) 



The time evolution of this vortex line is shown in Fig. 



VIII. VORTEX LINES IN RELATIVISTIC WAVE 
MECHANICS 

The description of vortex lines and their motion pre- 
sented in this paper is not restricted to nonrelativistic 
wave mechanics based on the Schrodinger equation. It 
can be extended to a relativistic theory based on the 
Klein-Gordon equation. Let us consider a relativistic 
wave function of a free particle in the form of a plane 
wave 



0k(r, t) — exp(ik-r) exp(— iu^t), 



(36) 



where — c-^/k 2 + (mc/h) 2 . This solution, in full anal- 
ogy with the nonrelativistic case, may serve as a gener- 
ating function for solutions with all kinds of vortex lines. 
For example, by taking first derivatives we generate a 
solution with one rectilinear vortex 



(x(t) cosx + iy(t) sinx) </>k(r, *)• 



(37) 



The only difference between this solution and the non- 
relativistic one is due to a different relation between 



the velocity and the wave vector. In the relativistic case 
it is v = Vk^k = Kk-I y/ (?ik/c) 2 + m 2 , instead of just 
Tik/m. Otherwise, the result is the same: the vortex line 
moves with the constant velocity v. The motion of two 
vortex lines can also be determined in the same way as 
in the nonrelativistic case. In the limit, when k — > 
the relativistic case reduces to the nonrelativistic case, 
Eq.©. 

With the help of the generating function we may also 
determine the motion of a vortex ring in the relativistic 
theory. The wave function that determines the motion of 
the vortex ring in the relativistic case (the counterpart 
of the nonrelativistic wave function (pl|)) is 



+ia(z(t) 



x(t) 2 + 
2h(l-c 2 (k 2 x 



y(t) 2 -R 2 



a,\/ (Kk/c) 



t) 0k(r,t). (38) 



The speed of the vortex depends on the arbitrary pa- 
rameter a and may exceed the speed of light when a is 
sufficiently small. This phenomenon, however, does not 
contradict the theory of relativity. It is quite similar to 
the so called superluminal propagation in optics p2| , p3[ , 
when the maximum of the signal travels with the velocity 
greater than c. In the present case it is the node of the 
probability density that may travel at an arbitrarily high 
speed. 

The presence of vortex lines is not restricted to rela- 
tivistic wave functions with only one component. They 
can also be found in wave functions with several complex 
components describing spinning particles. However, the 
freedom of attaching almost arbitrary vortex lines to any 
wave function is gone. There are stringent condition that 
must be satisfied in order to have physically admissible 
solutions. For example, a rectilinear vortex can be at- 
tached to a plane wave solution of the Dirac equation, 
but only when the direction of the vortex line is corre- 
lated with the particle momentum. 

Of particular interest are the vortex lines associated 
with photon wave functions. Simple states of that cat- 
egory are the well known multipolar states: the eigen- 
states of the angular momentum. The corresponding 
photon wave functions have rectilinear vortices running 
from — oo to oo, as in nonrelativistic wave mechanics. 
More elaborate vortex lines for photon wave functions 
are also possible but their study falls outside the scope 
of the present paper. 



IX. CONCLUSIONS 

In our study of the motion of vortex lines in wave me- 
chanics we have shown that this motion is determined by 
four elements: the shape of the vortex line, the shape of 
the wave function, the interaction between different vor- 
tex lines, and the external forces acting on the particle. 
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Complete isolation of these effects is not possible but one 
can gain some insight by appropriate limiting procedures. 
We found that quantum corrections to the vortex motion 
are responsible for the mutual interaction between vor- 
tex lines and also for the vortex creation and annihilation. 
We believe that the qualitative features of the vortex mo- 
tion will remain valid even in the much more complicated 
case of the mutually interacting particles, when the linear 
Schrodinger equation must be replaced by the nonlinear 
one. 
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FIG. 1. Creation of a vortex at a point. A vortex ring is created at a point where the wave function vanishes. 




FIG. 2. Instantaneous creation and instantaneous annihilation of a vortex pair with opposite circulation along a line. 
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FIG. 3. Two vortex lines evolving in time from negative to positive values of t. At t = these vortex lines form a pair of 
nonintersecting straight lines. 




FIG. 4. Precession of a rectilinear vortex in a magnetic field. 
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FIG. 5. Evolution of a vortex ring in a harmonic trap. 
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